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Abstract 

We prove that matrix Fredholm determinants related to multi time processes can be expressed 
in terms of determinants of integrable kernels a la Its-Izergin-Korepin-Slavnov (IIKS) and hence 
related to suitable Riemann-Hilbert problems, thus extending the known results for the single-time 
case. We focus on the Airy and Pearcey processes. As an example of applications we re-deduce a 
third order PDE, found by Adler and van Moerbeke, for the two-time Airy process. 
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1 Introduction and description of results 



Dyson, in [12], described how to implement a dynamics into random matrix models in such a way 
that the eigenvalues of the matrix behave like finitely many non-intersecting Brownian motions 
on the real line. In a suitable scaling regime we are lead to the study of certain interesting time- 
dependent determinantal random point processes generalizing the typical probability distributions 
appearing in random matrix models (see [21]). This paper is about two of these processes, described 
respectively by the matrix Airy and Pearcey kernels (see (2.2)-(2.4) and (3.2)-(3.4)). The Airy 
process was introduced in [19] and further developed in [15, 16], while the Pearcey process was 
introduced in [22, 18] in the context of non-intersecting Brownian motions and plane partitions. 

One important feature of the probability distributions appearing in random matrices it is that 
they are related to Fredholm determinants of integrable operators a la Its-Izergin-Korepin-Slavnov 
[14]. Namely, their kernels can be written in the form 

K( X ,y) = ^M (L1) 
x-y 

where /,<? are two vectors of a given dimension such that f T (x)g(x) — (so that, in particular, 
the kernel is non-singular on the diagonal). This simple fact has important consequences, as it 
reduces the computation of the Fredholm determinant to the study of a certain Riemann-Hilbert 
boundary-value problem canonically related to the operator's kernel (for a concise account see [13] 
and also the appendix A below). If, as it is very often the case, the determinant cannot be written in 
terms of elementary functions, still the Riemann-Hilbert formulation is very useful for finding some 
differential equations and studying asymptotic properties of the determinant. In this paper we will 
prove that the determinants of the Airy and the Pearcey matrix operators (namely the multi-time 
gap probabilities of the respective point processes) are equal to the determinants of some explicitly 
given integrable kernels (for details see Theorems 2.1 and 3.1) 

det(Id-x/A) = det(Id-i^A) (1-2) 
det (Id - XjP) = det(Id-i*Tp) (1.3) 

where A, P are respectively the matrix Airy and Pearcey operator, \i an d XJ are the indicator 
functions of a given collection of intervals and Ka, Kp are of the form 

K A (X,,) = ^^, ^(A,M)= fg( " )gp(M) . (1.4) 

A — fj, A — ^1 

Now f and g, for both cases, will be rectangular matrices of appropriate size, as in [13]. Our 
approach is the same used in [4] for the scalar Airy and Pearcey operators. We believe that this 
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method, with appropriate modifications, can lead to the proof that also others matrix kernels (like 
the Hermite, the sine and the Bessel described in [21]) are of integrable type. As an example of 
possible applications we describe how to obtain a system of isomonodromic Lax equations for the 
Airy and Pearcey processes. In the simplest case of the two-time Airy process on two semi-infinite 
intervals, we also re-derive from the Lax pair a partial differential equation originally discovered 
by Adler and van Moerbeke [1] and generalized for an arbitrary number of times by Dong Wang in 
[23]. 

Having a Riemann-Hilbert formulation for these Fredholm determinants can be used to address 
asymptotic of the multi-time Pearcey gap probabilities and how it connects with multi-time Airy 
gap probabilities, along the lines [4]. 

2 The Airy kernel 

The multi-time Airy process with times t\ < t% < . . . < t„ is governed by the matrix Fredholm 
operator A := A — B with kernels A(x, y), A(x, y), B(x, y) whose (i, j)-entries are given by [15] 4 

Aj(x,y) := Aij(x,y) - Bij(x,y) , 1 < i,j < n (2.1) 

J(x,ft)-6(y,\) 



(27Tl) 2 J 1r ^ J iR X + Tj-fi- T, 

9{x,n) := y-^- (2-3) 

1 (Tj-Tj) 3 (X-H) 2 (T-j-Tj )(»+«) 

Bii(x,y) := Xn< ——=— e 12 - - (2 . 4) 

V 47r ( T J - T i) 

Note that the matrix Bij(x, y) is strictly upper triangular. Here 7^ := jr — T{ and 7^ consists 
of two oriented rays, one from ooe^ to C and the other from C to ooe" 1 ^. C is a real number 
chosen so that C > maxj{rj}. Consider the collection of multi-intervals . . . ,I n } where 



T -S [4 ) ,af]U[af\af ) ]U...U[ap- 1 \af i) ] if % is even 
' < [af\af]U[af\af]U...U[ap\o ) if k 3 is odd. 



Moreover denote with \f its characteristic function. In the following, given any function h, we 
denote with h(dj) the column vector whose i th component is h(cij). Also we denote := 
diag(l, —1, 1, — 1, . . . , (— l) fe+1 ). The main result of the section is the following: 

Theorem 2.1 The following identity between Fredholm determinants holds 

det(Id - Xf A ) = det(Id - K A ) (2.6) 



4 Here the contours are slightly modified with respect to the one used in the reference. Still, our contours can be 
deformed into the original ones leaving the kernel of the operator unchanged. 
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where Ka is the integrable Fredholm operator acting on H := (@" =1 L 2 (iEL + Ti, C™)) © 1? (jr, C") 
with kernel 

fJ(A)gA(/i) 



AU(A,/i) = 



A // 



1 



f-A := :(/,A,l) • ■ • ' fA,n)t g : = ■ ■ ■ ,9A,n) 

s i e (0,A 3 ) X7jA) 

o fcl 



(2.7) 
(2.8) 



Am(A) := 



e SlAl x l R+r l (A) 



e- e(0 '^)x« +T ,(^) 
e e ^(o, ft )-«>, X7s(/1 ) 



(2.9) 



_ffere and below the variables A, /x, £ im£/i index, like for instance Aj, denote Aj := A — tj. 

Remark 2.1 XTie naming of Fredholm determinant in Thm. 2.1 is slightly abusive: strictly speaking 
here " det" is defined through the Fredholm expansion (B.l). In more precise terms the operator 
XfA is of trace-class, whereas XfB is Hilbert- Schmidt (this will appear clear in the proof below after 
(2.13)) and its kernel is patently diagonal- free. Therefore 



" det " (Id — A) := "det "(Id - A + B) = e -Tr A det 2 (Id - A + B) 



(2.10) 



which is the functional- analytically sound definition. Here det 2 denotes the Carleman regularized 
definition of determinant, for which we refer to Appendix B. 

Proof: We work on the entry of the kernel A(x,y) (2.1) and we observe that we can write 

(using Cauchy's theorem and after an appropriate shift of the variables of integration) 

A i j(x,y)xi i {x) = 



d £.^(_)*+i e &(a 4 w -s) x 



2iri 



dX 



dfi e e ( a i ) .w)-fl(o.A j )+»A J 

2ni (C-^(M-A) + Xt ^ 



dfi e e( - ' >t J -i)- e ( 'N)+VN 



2iri 
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We deduce that XfA = T x &aT where T = diag(7i, . . . , T n ) is the following diagonal Fourier 
transform, \ij := fj, — Tj 



Tj : L 2 (iR + Tj) 



L 2 (R) 
1 

V2hr 



h(x) i y 



L 2 (iR + Tj) 
1 



'2m 



< x h(x)da 



and &a has kernel explicitly given by entries 



(£aM£,a) = ]T(-i/ +1 



1=1 



2m (t-Mn-X) 



+ Xr i <Tf 



.(2.11) 



Here &a is an operator 5 acting on L 2 (iM. + Tj, C n ). Using (2.11) we can write 

&A = Qa ° -Fa + Ha where Q a-, Fa an d Ha are the matrix Fredholm operators 



J A :0L 2 (iR + r„( 

n 

H A :($L 2 {iR + T t , 



L 2 h R ,C n ) 

n 



g A :L 2 ( lR ,C n )^®L 2 (i 



i=l 



with kernels having entries given respectively by 

e |e(o, w )-0(o,A 3 ) 



(JU.)»i(M,A) := 



/i — A 

(«A)««,A) ^x^^-l)^ 



(aA)««,M) = *«E(- 1 ) <+x 



£=1 



1 fl(a^ ) ,A < )-e(0,Af)+a^ ) f 4 



e-A 



(2.12) 



Now consider the Hilbert space i? := (0" =1 L 2 {iR + t u C™)) © i 2 (7^, C") and : H — > H the 
matrix Fredholm operator written in the matrix form (naturally-related to the splitting of H into 
its two main addenda) 



K. 








Qa 





(2.13) 



It is easily verified by the convergence of the L 2 norm of their respective kernels that all the operators 
J~a,GatHa are Hilbert-Schmidt and thus XfA — 1~~ 1 Ga J~a7~ is of trace-class being unitarily 
equivalent to the product of two HS operators. On the other hand and by the same argument 
XfB = —T~ 1 'HaT is provably only Hilbert-Schmidt. 

5 It is worth noticing that as written in (2.11), is already an integrable kernel a la IIKS. Still, in order to have 
a more "well behaved" Riemann-Hilbert problem, we go on manipulating the kernel a little bit more. 
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To avoid confusion we shall temporarily denote the determinant defined by the Fredholm ex- 
pansion (B.l) by "det". Then "det"(Id — K a) = det2(Id — Ka) since Ta,Ga and Ha ah have 
diagonal- free kernels. 

Let's also define an other operator K' A written as 



K' A 















(2.14) 



This operator K' A is only Hilbert-Schmidt (not trace-class) and hence some care needs to be paid 
when inserting it into a (regularized) Fredholm determinant. To this end we start observing that 
its Carleman determinant det2 (see App. B) is well defined and det2(Id — K' A ) = 1. Then we have 
the following equalities (the operator A is defined in (2.1)): 

"det"(Id - K A ) = det 2 (Id - K A ) = det 2 (Id - ^)det 2 (Id - K' A ) = 
= det 2 (Id - ^)e- Tr(e ^°^ ) = dct 2 (Id - xrA)e~ TrA = "det "(Id - X fA). ( 2 - 15 ) 

The last equality in the first line is just an application of (B.3); then we go from the first line to 
the second using the unitary Fourier transform T and the last equality is due to the fact that the 
kernel of the operator A — A is diagonal-free (see Remark 2.1). It just remains to prove that Ka 
has kernel defined by the equations (2.7)-(2.9). In view of (2.13) and definitions of Fa, Ga, Ha this 
is just a matter of straightforward computations. Q.E.D 

2.1 Riemann— Hilbert problem for the Airy process 

Here we show how, using Theorem 2.1, we can relate the computation of the Fredholm determinant 
of the matrix Airy operator to the theory of isomonodromic equations. We start defining the 
following RH problem that, according to the IIKS theory, is naturally related to the Fredholm 
determinant det(Id — Ka)- 

Problem 2.1 Let T(A) be the sectionally analytic function (if existing) that solves the RHP 

n 

T+(A) =r_(A)(l-27rzG A (A)) , Xe 1 :=^ R Li\J(tR + T l ) (2.16) 

i=l 
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-6>(5i,Ai) 



„— 9(a2,A2) v , 
e A! 



C 2 ,i(A) 







Ci>j (A) = x„+r 4 (A) f ( -i)*+i e - e («r,A 4 ) + e(ar.^)) 



C n ,i(A) 



€ Matfc x kj) 
=i...fci.t=i...fc,- 



Remark 2.2 TTie jump on the contours i$L + Tj is oscillatory on the first column due to the term 
e -0(o 3 -,A 3 -) while on the other columns is exponentially decaying (this is due to the fact that the leading 
term of the exponential in Cj j is 3(Tj — Tj)\ 2 with i < j). In order to have exponential decay in 
all columns we can now deform all the contours iM. + Tj to jl + Tj with := — 7^. We could also 
have collected all these new jumps on the single contours by multiplication of the jumps. The 
resulting jump, however, would be slightly more complicated. 

It is straightforward that Gyi(A) = f,4(A)gT(A) and that g^(A)fU(A) = 0„. Hence, using Theo- 
rem A.l, we arrive to the following 

Theorem 2.2 The Fredholm determinant det(Id— XT 'A) is equal to the isomonodromic tau function 
tjmu W (A. 6) related to the RH problem (2.1). In particular Vz = 1, . . . , n and = 1, . . . , fej we 
have 



dA 



aindet(Id-xM) = J Tr(r: 1 (A)r'_(A)S a (A)^ 

3 a (A) := (9(1 - 2ni G A ) (1 - 2m GaY 1 = ~2m dG A (l + 2m G A ) 

where we have denoted by d any of the derivatives 3 w , d Ti . 

Theorem 2.2 implies also some more explicit differential identities by using the Miwa-Jimbo- 
Ueno residue formula; note first that the jump matrices G a can be written as 



G A (X) 



T A (\)Q(0) e -T A (X) 



(2.19) 



(2- 



where G? is a constant matrix (consisting of only ±1 and 0) and 
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T A (X) - diag (if, f . . . , f JO) , Tf := * £ £ 9(af\X l ) (2.20) 

ff := T^l - diag (9(a^,Xi), . . . , 6{af l) , A,)) , Tr T A (A) = 0. (2.21) 

Therefore the matrix ^ a{X) := r A (A) c Ta ^ solves a RHP with constant jumps and hence is 
(sectionally) a solution to a polynomial ODE 

9 A * A (A) = L(A)* A (A) , degL(A) - 2. (2.22) 

It was shown in [3] and it is not hard to see using Cauchy's residue theorem (see also [5]) that 

/ Tr(ri 1 (A)r'_(A)S a (A)) |4 = - res Tr (T- 1 (X)T'(X)dT A ) dX (2.23) 

where the (formal) residue in (2.23) simply stands for minus the coefficient of the power A -1 in 
the asymptotic expansion of the argument. Direct application of (2.23) using the expression of Ta 
(2.20) and Theorem 2.2 yields the following proposition. 

Proposition 2.1 The Fredholm determinant det(Id — Xf ~A) satisfies 

9 a wlndet(Id-xp4) = -( r i)i+(£ 3<s fc j )+^,i+(E j<i k s )+i (2-24) 

ki 

d Ti Indet(Id-xM) = £(2r i r 1 +T\ - 2r 2 ) 1+(E . <4 ^ )W+(E . <s (2-25) 
l=i 

where Pi,r2 and T(X) are the coefficients of the asymptotic expansion of the solution of the RHP 
2.1. 

Remark 2.3 (On the existence of the solution to Prob. 2.1) According to Thm.2.1 and fol- 
lowing the general IIKS theory of integrable operators, it follows that the solution to Problem 2.1 
exists if and only if the Fredholm determinant det(Id — Ka) does not vanish. The latter determinant 
is equal to det(Id — XfA) by Theorem 2.2, and this follows from the probabilistic interpretation of the 
same. Of course this should be proved also independently, for example by showing that the operator 
A has operator norm bounded by one. 

2.2 Lax formulation of a PDE for the Airy process 

Using the representation (2.22) above we can give an independent verification for the nonlinear 
PDE satisfied by G(r,a,b) = logdet(Id — X{(a.oo).(fc,oo)}^4) as found by Adler-van Moerbeke in [1] 
and generalized by D. Wang [23] to the case of an arbitrary number of times (here r = t-i — t\): 

{^d w - Wd E ^j (d% -d 2 w )G + 2rd 3 TEW G = {d 2 EW G, d 2 E G} E (2.26) 
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where E = 9 4^, W = and {f,g}E '■= d E f g — fd E g- We briefly describe how to use our 
setup to verify the equation (2.26); the details of the computations involve a significant amount of 
completely straightforward algebra and will not be reported here. We are in the case n = 2 and 
ki — &2 = 1 so that the RHP is 3 x 3; since the process is stationary we can assume t\ = and 
T2 = r > 0. The matrix ^ A (X) — r a(^)g Ta ^ solves a RHP with constant jumps, where 

(a, A) + 0(6, A - t) 9(b, A - t) - 20(a, A) 0(a, A) - 20(6, A - r) ' 



T A (A) = diag 



(2.27) 



Denoting by the coefficient matrices in the expansion of the solution as in Problem 2.1, one can 
write a set of compatible PDEs 



d w <£ A {\)=W(\)* A (\) 



d x y A (\) = L(\)y A (\) 



(2.28) 

d T * A (\) = T(A)* A (A) (2.29) 



with all the matrices polynomials in A of degrees 2,1,1,2 re- 
spectively. From the equation (2.28) one can express all the co- 
efficients Tj in terms of the coefficients of Ti and the entries 
(1, 2), (2, 1), (3, 1), (1, 3) of T 2 (we shall refer to these coefficients 
as the "generators"). On the other hand, equating the expansions 
of equations (2.29) allows to express all derivatives in terms of these 
generators uniquely. From the specialization of Prop. 2.1 we obtain 

d E G = - Tr(r!diag(0, 1, 1)) , 3 W G = - Tr^diag^, 1,-1)) 
d T G = Tr ((2tTi + Tj - 2r 2 )diag(0, 0, 1)) 

Then all the higher partial derivatives can be straightforwardly ex- 
pressed in terms of polynomials in the generators. By simply plug- 
ging these lengthy expressions into the equation (2.26) yields an 
identity. 




Figure 1 : The contours for the 
Pearcey kernel and Riemann- 
Hilbert problem 



Remark 2.4 Although this method gives a straightforward verification of the equation, it is not 
constructive in that it does not provide an effective method of finding an equation. In general we 
would expect as many nonlinear PDEs as the number of independent variables. 



3 The Pearcey kernel 

In this section, slightly modifying the approach used for the Airy kernel, we show that also the 
determinant of the matrix (multi-time) Pearcey operator can be expressed as a determinant of a 
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related integrable operator. The matrix Pcarccy kernel (sec for example [22]) has (i, j)-entry given 

by 



Qi]{x,y) 



■= Pij(x,y) - Q(x,y 



' 3 







d/i / dA 



(2ttz) 2 

M 4 



>i 2 



4 2 

Xr,< 



/i — Xfi. 
1 



T ' / ' 

3 V 27r ( T i _ T i) 



"2(t,-— r;) _ 



Xn<7 



dX 
2ri 



(3.1) 
(3.2) 

(3.3) 
(3.4) 



where the contours 7l,7p are indicated in Fig. 1. Here we consider the Fredholm determinant 
det(Id— XfP) 6 where Xf' is the characteristic function of the collection of multi-intervals { J\, . . . , J„} 
with Jj :— [a"p ,a^] U [aj 3 \aj 4 ^ U . . . U [aj 2fcj_1 \ a^ 2fe "] (note that for the Pearcey case we do not 
have semi-infinite intervals). The analog of Theorem 2.1 is the following. 

Theorem 3.1 The following identity between Fredholm determinants holds 

det(Id - XjP) = det(Id - K P ) (3.5) 
where Kp is the integrable Fredholm operator acting on V := L 2 (iM. U U jl) with kernel 
f£(A)gp(M) 



A — n 

(/p,1j ■ • ■ ; /p,ra)i 



(SP,1> ■ ■ ■ i9P,n) 



(3.6) 
(3.7) 



(/p)i(A) := 



e^ (0 < A) X 7 ,u 7 ,(A) 



e SlA X 4 R(A) 



J2fe„ 



£2fcie 



-Si^+(Tj-ri) i 



"XiR^) 



-aj_i^+(Tj-rj_i)' 



"XiR(/") 



3 -e^ (0 ^- s ^X 7i u 7R (/i) 



J2fc„ 



(3.8) 



3 Note that we have det(Id - X /P) = dct(Id - x jPx f) 
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Remark 3.1 The same caveats about the correct definition of det as in Remark 2.1 apply here: 
as a matter of fact, a small modification to the kernel would turn it into a fully trace-class opera- 
tor (without changing the expansion of its "Fredholm" determinant), but for reasons of space and 
uniformity with the previous treatment we use the same approach. 



Proof: 



We work on the (i, j)-entry of the kernel and we observe that we can write 

2k, 



p ij(x,y)xji{x) 



dX 



du e e *( a > f) ^)-e 3 (o,A)+yA 



27r * J 1L u 1H 2 ™ (£-m)(m-a) 



dfx e 



fa— Ti)4— 0} 'fl + yfJ, 



2ni 



(3.9) 



Then 



X,fP = T^KpT 
where T is the Fourier transform (acting diagonally) 

T:L 2 (iM) — -> L 2 (R) r _1 :L a (R) 
1 



L 2 (iR) 



e^/(/x)d/x 



/2l7T 



e~ flx h(x)dx. 



(3.10) 



and ftp has kernel explicitly given by entries 



2k, 



(*p)«(£,A) = £(-i) m 



^ e e i (af ) , M )-e i (0,A)+af ) £ e ( T;f _ n) ^._ W (A _^ 

tlu 7 « 27r * Xti<Tj £-A 



(3.11) 



.ftp is an operator 7 acting on L 2 (?]R, C n ) and such that .ftp ~ Qp o Tp +Hp with 

Jp : L 2 («, C") — > L 2 (( 7L U 7 p, C n )) , £ P : L 2 (( 7l U 7fl ), C n ) — ► L 2 («, C") 
Hp : L 2 (iR,C n ) — ► L 2 (iR,C n ) 

and (£, j)-entries given respectively by 



2k, 



(J»y(^,A). 

/i — A 

(«p)«(£,A) :=Xn<T .^(-iy 



(^) ii ^A*) = *«D- 1 ) /+1 



e ar ; (C-A) + (r J -r l )4 r 



C-A 



(3.12) 



7 It is worth noticing that £p as written in (3.11) is already an integrable kernel a la IIKS. In particular and most 
importantly, there is no singularity on the diagonal f = A for the second term thanks to the alternating sum. 
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Now consider the Hilbert space V := L 2 («,C n ) L 2 (( 7L U 7 R ),C n ) and K P : V — > V the 



matrix Fredholm operator written in the matrix form (naturally-related to the splitting of V) 



Kp 








Qp 


Up 



(3.13) 



Here again, as in the case of the Airy operator, all the operators FpjGpjUp are Hilbert-Schmidt 
so that det2(Id — Kp) is well defined. The very same argument used for the Airy operator leads 
us to conclude that det(Id — XfP) — det(Id — Kp) (we recall once more that here "det" stands 
for the Fredholm expansion (B.l)). Formulas (3.6)-(3.8) are again a matter of straightforward 
computations using definitions of Tp,Qp,U.p and (3.13). Q.E.D 

3.1 Pearcey process and isomonodromic equations 

In this section we relate the computation of the determinant of the matrix Pearcey operator to the 
theory of isomonodromic tau functions, thus potentially leading to a Lax representation. We start 
defining our RH problem 

Problem 3.1 Let T(A) be the sectionally analytic function (if existing) that solves the RHP 
r+(A) =r_(A)(l-27TzG P (A)) , A e 7 := lR U 1l U *R 



r(A) = 




Gp(A) := 



-e 2 (3 2 ,A)„ 



A — > oo 



D 2) i(A) 



(3.14) 
(3.15) 





D i,i(^) = X«(A) 



■ e " (3 "' A) XiR D„ 4 (A) 

t+l p (a< s) -a< t) )A+^A 2 



(-l) t+1 e 



■ • ■ -Dn,n— 1 

e Mat(2fej x 2kj) 

= 1...2ki,t=1...2ki 



(3.16) 



Remark 3.2 In contrast to the Airy case here the jumps of the RH problem 3.1 are already expo- 
nentially decaying and no deformations are needed. 

Remark 3.3 The same considerations as in Remark 2.3 about the existence of the solutions for 
the RH problem 3.1 apply here. 
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The reader can verify that Gp(A) = fp(A)gp(A) and that fJ(A)gp(A) = 0„ (see (3.8)) so that, 
similarly as in the previous section, we obtain the following theorem. 

Theorem 3.2 The Fredholm determinant dct(Id — XfP) equals to the tau function tjmu i n (A. 6) 
related to the RH problem (3.1). In particular Vz — 1, . . . , n and W £ = 1, . . . , 2ki we have 

dlndet(Id- Xr P) =U(x- 2 rtG P )(d) := / Tr(r: 1 (A)r'_(A)S a (A)) ^ 
S 9 (A) := (5(1 - 2m G P )) (1 - 2m Gp)" 1 = -2m <9G P (1 + 2m G P ) 
where d stands for any of the derivatives d (*) , d Te . 

Again using the Miwa-Jimbo-Ueno residue formula we get some more explicit differential equations 
in terms of the coefficients of the symptotic expansion at oo. As in the Airy case the jump matrix 
Gp can be written as Gp(A) = e Tp ' A 'Gp'e" T ' p ^ where G^ is a constant matrix (consisting of 
only ±1 and 0) and 

1 n 2ki 

Tp(A) = diag (T p °\fP, . . . ,2*0) , TP ;= \ £ £ 6,(af\ A) (3.17) 

f« :=T| 0) l-diag(e j ( a f ) ,A),...,e(af' £ ' ) ,A)) , Tr T P (X) = 0. (3.18) 

The matrix 4'p(A) :— T(X)e Tp ^ solves a RHP with constant jumps and hence is (sectionally) a 
solution to a polynomial ODE. As in the Airy case the integral in the Theorem 3.2 can be converted 
in a formal residue, namely 

J Tr(r: 1 (A)r'_(A)S (A)) ^ = - re^ Tr (T- 1 (X)V (X)dT P ) (3.19) 
We then find exactly as for Prop. 2.1 

Proposition 3.1 The Fredholm determinant det(Id — XfP) satisfies 

lndet(Id-XfP) = -(^lh+^^i^+iA+CE^i^+e ( 3 -20) 

^ 2ki 

5 Ti lndet(Id- X/ P) = -£(r? -2r 2 ) 1+(E . <i2fc . )+ , il+(E . <i2fc . )+ , (3.21) 

where ri,r2 and T(A) are given by the solution of the RHP 3.1. 
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A Integrable kernels and isomonodromic tau functions 



In this appendix we review some basic facts about integrable kernels [14] and, for sake of com- 
pleteness, a theorem we originally proved in [4] (see also [5]). Given a piecewise smooth oriented 
curve C on the complex plane (possibly extending to infinity) and two matrix-valued functions 
f , g : C — > Matpxfe(C) we define the kernel K as 

f T (A)g( M ) 
K{X,n) := — . 

A — [1 

We say that such kernel is integrable if f T (A)g(A) = 0^ (so that it is non-singular). We are 
interested in the matrix-valued operator K : L 2 (C,C k ) —> L 2 (C,C k ) acting on A:-vector functions 
via the formula 

(Kh)(X) = J K(\,fj,)h(fi)d[j, 

and, in particular, we are interested in the Fredholm determinant det(l — K) defined as in (B.l). 
Denoting with d the differentiation with respect to any auxiliary parameter on which K may depend, 
we have Jacobi's formula 

d\ogdet{l-K) = -Tr((ld + R)dK) (A.l) 

where R is the resolvent operator, defined as R = (1 — K)^ 1 K . Moreover R is again an integrable 
operator 

F T (A)G( M ) 
R{\n = ; 

A — f-l 

where F, G are given by F(A) := T(A)f(A), G(A) := (r- 1 ) T (A)g(A), and T is the solution of the 
following RH problem: 

r+(A)=r_(A)M(A) XeC (A.2) 
r(A)~ l + C(A- 1 ) A^oo (A.3) 
M{X)= 1 - 2nif (A)g T (A) (A.4) 

The theory also guarantees that the solution to the above problem exists if and only if the Fredholm 
determinant does not vanish. Now suppose that the kernel the matrices f , g and thus K (and hence 
the Riemann-Hilbert problem (A.2)-(A.4)) depend smoothly on parameters. On the space of these 
deformation parameters, we introduce the following one-form 8 

dA 



u M {d)~ / Tr(ri 1 (A)r'_(A)S a (A))^ , E d (X):^dM(X)M- 1 (X). (A.5) 



8 Here and below we will denote with ' the derivative with respect to A. 
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Moreover, in the case that ujm (which is a condition only on the jump matrices as seen in [3] is closed 
in the space of the deformation parameters, it is also defined (at least locally) the isomonodromic 
tau function tjmu such that 



d log t j mu ■= w M (<9) 



(A.6) 



for every vector field d in the space of deformation parameters. The definitions (A. 5), (A. 6) are 
posed for arbitrary jump matrices; in the case of the RHP (A.2)-(A.4) tjmu in (A.6) and the 
Fredholm determinant are related as in the theorem below. 

Theorem A.l ([4]) 9 Let f (A; s), g(A; s) :CxS — ► Mat pX fe(C) and consider the RHP with jumps 
as in (A.2)-(A.4). Given any vector field d in the space of the parameters S of the integrable kernel 
we have the equality 



where u>m{9) is as in (A. 5). 

In the cases we treat in this article the additional terms in (A.7) vanish because of the particular 
shape of the matrices f , g: indeed, due to their particular shape, a stronger condition holds that 
f T (A)g(^) = Ofc when A and /i both belong to the same contour. Hence the tau function tjmu, 
defined in (A.6), coincides with the Fredholm determinant det(l — K). 

B A brief reminder of regularized Fredholm determinants 

We refer to [20] for the relevant details: we shall need only the elementary facts which we recall 
here. In general the Fredholm determinant of an operator of the form (Id — G) can be defined 
only when G is of trace class. If G is represented as an integral operator on a (separable) Hilbert 
L 2 {X, d/i) with kernel G(£,r)) (we abuse notation by using the same symbol for the operator and 
kernel) then 10 



The trace ideals I p , p £ N consist of operators whose p-th. power is trace-class [20]; in particular I2 
consists of Hilbert-Schmidt operators. For G £ I p one can define following Carleman a regularized 

9 Actually the article [4] treats the case k = 1, but the proof does not change considering this more general case. 
10 Note that, if X = {1, . . . , n} X S and fi = v % A where v is the counting measure, the formula above covers the 
case of matrix-valued kernels on S equipped with measure A, as in [15]. 



u M {d) = 91ndet(l-if)+ I <9Tr(f ,T g)dA + 2iri [ Tr(g T f'<9g T f)dA 



(A.7) 




det(Id - G) := 1 + Y, / det[G(6, &)] w < fc J] d M (&) 




(B.l) 
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determinant det p (Id — G) which has the same main property of vanishing iff the operator is not 
invertible. In particular for Hilbert-Schmidt operators one has 

dct 2 (Id-G):=l + V- / det[Gte,0)(l-^)] i , j < fe Tl d ' u ^) ( B - 2 ) 
fc=i k ' Jxk i=i 

that is, we omit the diagonal elements in the determinants under the integral sign. It satisfies 

• if G e Ii n 2 2 then det 2 (Id - G) = det(Id - G)e TlG 

• if Gi,G 2 are Hilbert-Schmidt operators (and hence GiG 2 is trace class) then 

dct 2 (Id - Gi)det 2 (Id - G 2 ) = dct 2 (Id - Gi - G 2 + GiG 2 )c Tr ( GlG2 ) (B.3) 

An interesting occurrence (which is used in this article) is that if G is just HS but its kernel 
vanishes on the diagonal G(£, £) = then the series defining det 2 (Id — G) is identical to the regular 
dct(Id — G). The reason for still wanting to distinguish det 2 from det in this case is simply that G 
may fail to have a trace and in a different basis of the Hilbert space the ordinary det may simply 
be ill-defined. 
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